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We obtain an exact formula for the equilibrium free energy of a charged quantum particle moving 
in a harmonic potential in the presence of a uniform external magnetic field and linearly coupled 
to a heat bath through the momentum variables. The heat bath is modelled as a collection of 
independent quantum harmonic oscillators. The derived free energy has important differences from 
that in the case of coordinate-coordinate coupling between the particle and the heat-bath oscillators. 
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I. INTRODUCTION 

Consider the dissipative dynamics of a charged quan- 
tum particle moving in a harmonic potential in the pres- 
ence of a magnetic field and coupled to the external en- 
vironment. Such a situation is frequently encountered in 
typical experiments designed to study the magnetic re- 
sponse of a charged quantum particle in the context of, 
e.g., Landau diamagnetism jlij, quantum Hall effect Q], 
and high-temperature superconductivity [3| . 

In a recent paper [J| , we analyzed the dynamics of the 
charged particle by regarding the environment as a quan- 
tum mechanical heat bath or reservoir, and by invoking a 
gauge-invariant system-plus-reservoir model. We consid- 
ered the heat bath to consist of a collection of indepen- 
dent quantum harmonic oscillators, while its interaction 
with the charged particle was modelled in terms of bi- 
linear coupling between the momentum variables of the 
particle and the oscillators. The magnetic field was taken 
to be uniform in space. In that work, we obtained the 
quantum Langevin equation (QLE) for the reduced dy- 
namics of the particle, in which coupling to the bath is 
described solely by an operator- valued random force and 
a mean force characterized by a memory function. The 
magnetic field appears in the QLE through a quantum- 
generalized classical Lorentz force term, while the ran- 
dom force is independent of the field. These aspects 
are also present in the QLE for the case of particle-bath 
coordinate-coordinate coupling [5|. However, there are 
important differences, e.g., in the form of the memory 
function which now depends on the magnetic field, in the 
expressions for the symmetric correlation and unequal- 
time commutator of the random force, and also in the 
fact that the inertial term and the harmonic potential 
term in the QLE get renormalized by the coupling con- 
stants. 

In this sequel paper, we further our study of the 
same system by deriving an exact formula for the free 
energy of the charged particle in thermal equilibrium. 
Our result shows important differences in the form of 
the free energy with respect to that for the coordinate- 



coordinate coupling @- Thus, this work, in addition to 
the previous one [J| , underlines the subtle role played by 
the form of the coupling (coordinate-coordinate vis-a-vis 
momentum-momentum) in characterizing the dynamics 
and thermodynamics of the charged particle. 

The paper is structured as follows. In the following 
Section, we define the system of study and write down 
the equations of motion for the charged particle and the 
heat-bath oscillators. Section Hill forms the main part of 
this paper in which we present a detailed derivation of the 
equilibrium free energy of the charged particle, and point 
out its differences from the case of coordinate-coordinate 
coupling. 



II. SYSTEM OF STUDY AND EQUATIONS OF 
MOTION 

The system considered in Ref. ;4i] comprises a charged 
quantum particle moving in a harmonic potential in the 
presence of a uniform external magnetic field, and lin- 
early coupled to a quantum heat bath through the mo- 
mentum variables. The heat bath is modelled as a collec- 
tion of a large number of independent quantum harmonic 
oscillators. The Hamiltonian of the system is 
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where e, to, p, r are, respectively, the charge, the mass, 
the momentum operator and the coordinate operator of 
the particle, while wq is the frequency characterizing its 
motion in the harmonic potential. The jth heat-bath os- 
cillator has mass mj, frequency Uj, coordinate operator 
Qj, and momentum operator pj. The dimensionless pa- 
rameter Qj describes the coupling between the particle 
and the jth oscillator. The speed of light in vacuum is 
denoted by c. The vector potential A = A(r) is related 



to the uniform external magnetic field B = {Bx,By,Bz) 
through B = V X A(r). The field has the magnitude 

B=,jBl + B^y+Bl 

The commutation relations for the different coordinate 
and momentum operators are 



[r^Pp] = iflSaP, [qja,Pkl3] = iflSjkSa/S, 



(2) 



while all other commutators vanish. In the above equa- 
tion, Sjk denotes the Kronecker S function. Through- 
out this paper, Greek indices (a, /3, . . .) refer to the three 
spatial directions, while Roman indices (i, j, fc, . . .) repre- 
sent the heat-bath oscillators. Moreover, we use Einstein 
summation convention for the Greek indices. 

The Heisenberg equations of motion for the charged 
particle are 
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where m,- is the "renormalized mass" , 
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and dots denote differentiation with respect to time. The 
Heisenberg equations for the heat-bath oscillators are 
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These equations lead to 
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For details on the derivation of Eqs. (|31) - ®> see 



III. EQUILIBRIUM FREE ENERGY 

In order to obtain the equilibrium free energy of the 
charged particle, we start with computing the mean en- 
ergy {Hq) of the system ([T]) in thermal equilibrium at 
temperature T. Using Eqs. ^ and (O, we get 
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Thus, to find (Hq), we need equilibrium averages, such 
as (r^) and (r^), which may be obtained by knowing the 
equilibrium autocorrelation function of the position of 
the charged particle, defined as 

r;At-t')^l{r,it)rAt') + rAt')r,it)). (11) 

Similarly, (q^), (Qj), and (q^ • q^) are obtained from the 
position autocorrelation function of the heat-bath oscil- 
lators, defined as 



V'"? a -t') = - 



{qjp{t)qkAt') + qkAt')qjp{t)), {12) 



while (r-qj-Hqj-r) may be computed from the correlations 
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Cp.(i - 1') ^ o{rp(.t)qMt') + qMt'Vpit)): (13) 



^TpAt ~ t') EE -{qj^pity^t') + r,(i')g,,p(0)- (14) 



Correlations such as in Eqs. pT|) - (|14p are computed 
conveniently by invoking the fluctuation-dissipation theo- 
rem that relates such equilibrium correlations to response 
of the system to small external perturbations . To this 
end, consider a weak external force f (i) to act on the 
charged particle and another set of weak external forces 
{fj{t),j = 1,2,...,A^} to act on the oscillators. Here, 
f(t) and ij{t) are c-number functions of time. We take 
the perturbed Hamiltonian to be of the form 
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In the presence of these external forces, the Heisenberg 
equations of motion ([3]) and © remain the same, while 
those given by Eqs. ^ and ([7]), respectively, are modified 
to 
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Consequently, Eqs. (|5]) and ©, respectively, are modified 
to 
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r = -mujgr+-{rxB)+2_^gjm,ujjqj-2^ -fj+f, Using Eq. (gH) in either Eq. (gHl) or Eq. dm gives 
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mjC{j = —'mjUj'jqj+gjmLL!Qr~'-^{rx'B) + tj—gj{. (19) 

Takiiig the Fourier transform of Eqs. ([T^ and ([TO|. 
we get 
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where epari is the Levi-Civita symbol. 
Equations ^ and dH]) give 



lujgjt 



^paii-^i1 I ^1 



(21) 



N 2 

^P.(^)^-. = G(^)/p + J2 zrf?^/... (22) 



^mj(w^. -u;2) 



where 



AT 



G(a.) = l-5:.^^^-^''"^"^' 



,■=1 J^ J ' 



(23) 
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and 
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From Eq. ([22]), we get 
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The functions Q:p-y(a;), /3j^p-y((u;), Aj^p^(a;), and jjk.p-yit^) 
are coefficient matrices of the response of the system to 
the external perturbations, f(i) and {j(t), and may be 
interpreted as generalized susceptibilities. 

We may now use the fluctuation-dissipation theo- 
rem which relates the Fourier transform '0pj(aj) of 
the position autocorrelation function i^'Zrit — t') with 
G{(jj)apa{ijj): 



Vpl{^) = ^cot\i\^^]G{u:)[ap„{uj + iQ+) 



2i 



'-up 



2kBT 



(33) 



where * denotes complex conjugation (Note that G{uj) is 
real). Similarly, one has 



^llpA-) 



2i 



■ coth 



r huj 



2kBT 



Al]k,pa{^ + iO+) - l*k,ap{^ + *0+)], 



Cp-(^) = ^C°tl^ 



huj 



and 



Cp-(^) = ^c°t^ 



2kBT 



hcu 



2kBT 
In passing, we note that 
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which imphes that Im[acrp(w)] is an odd function of w, 
while Ile[a^p{uj)] is an even function of w. These proper- 
ties are also shared by I3j^a-p{^), ^j.api^), and jjk.api^)- 
In the following, we use Eqs. (|33)) - (p6| to evaluate 
the various averages involved in computing (-ffo)- 



A. Computation of various averages 

1. Averages involving r 



Differentiating Eq. (|4T|) successively with respect to t 
and t' , and finally putting t — t' , we get 
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We start with decomposing ap^ into symmetric and 
antisymmetric parts as 



ap^iu) = al„{uj) + al„{uj). 



We then have 
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= [a;,{u:) - a^pA^)*] + [a^,(w) + a^^iuj)*] 

= 2tlm[a%iw)]+2Re[a;,{uj)]. (39) 

Using Eqs. dTT]), dM]), ^, and the properties that 
Im[apcr(w)] is an odd function of w, while both G{lu) and 
Re[ap^(aj)] are even functions of w, we get 



duj G{iu)lin[a%iu + iO+)] 



X coth 



huj 



cos[a;(i - t')] 



.2kBT\ 

duj G'(a;)Re[a^^(w + iO+)] 
r tujj 



X coth 



2kBT 



sin[a;(i-i')], 



(40) 



which for p ~ a gives 
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Then, putting t = t' , we get 
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where from Eq. ()27|). we have 
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One may proceed similarly to the preceding subsection 
to obtain 
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where ^jk.pp{uj) is obtained from Eq. ([32]) as 
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We have 
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Adding Eq. (|5T|) to Eq. dSO]), we get 
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Differentiating Eq. (|5^ successively with respect to t 
and i', and then putting t = t' , we get 
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Using Eqs. ([35]) and ([36]) for V'j,pcr('^) and i^f^^ p{ijj) , re- 
spectively, and noting that the imaginary part of both 
Pj.paiy) and Aj p£r(aj) is an odd function of w, while their (ITUl) . we get 



B. Mean energy 
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The mean energy of the heat bath is 
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Here, in obtaining Eq. ([57l) . we have used Eq. (j46| . w + iO+, and by using the result that 
To arrive at Eq. (|58p . we perform the first integral on 



the right hand side of Eq. (|57l) by noting that to in the 
integral is approached from above the real axis, w 
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with P denoting the principal vahie. Note that the first We define the particle energy Uo{T,B) as the mean 

term on the right hand side of Eq. (1551) is the mean energy of the system of the charged particle interacting 

energy Ub{T) of the heat bath in the absence of any with the heat bath minus the mean energy of the heat 

coupling with the charged particle: bath in the absence of coupling with the particle: 
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where we have used Eqs ((29)) . (|46l) . and (1551) . After simplification, we get 
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From Eqs. ^ and ([Ml), we have 
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From Eq. ([25]) . we get 
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Now, using Eqs. (|M1), dSS]) and §^ in Eq. dHS]), we 
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where u{u!,T) is the Planck energy of a free oscillator of 
frequency lo: 
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C. EVee energy 

The free energy Fq (T, B) of the charged particle is ob- 
tained from Eq. ([S7)) . since the two quantities are related 
as Uo{T, B) = Fo{T, B) - T^^MI^. We get 

Fo{T,B) = - f duj f{uj,T) 

xlml ^ ln[Det a{uj + iO+)] 

iDet a(w + iO+) 



.A( 



duj 
ujeB\^/dG^ 
V du) 



(69) 



where f{uj,T) is the free energy of a free oscillator of 
frequency w: 
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Equation (|69|) is the central result of the paper. To 
make explicit the contribution of the external magnetic 
field to the free energy, we write Eq. (l63t as 
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is the susceptibility in the absence of the magnetic field. 
Using Eq. dTT]) in Eq. i^, we get 
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is the free energy of the charged particle in the absence 
of the magnetic field. The contribution from the latter is 
contained in the two terms AiFo(T, B) and A2Fo(T, B), 
given by 
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Let us now comment on the form of the free energy 
([75]) with respect to that for coordinate-coordinate cou- 
pling between theparticle and the heat-bath oscillators, 
obtained in Ref. @ ■ In the latter case, the free energy is 
given by 



Fo(T,B)==Fo(r,0) + AiFo(r,B), (77) 



where Fo(T, 0) has the same form as in Eq. ([7^ . while 
AiFoiT,B) is given by 
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We thus see that with respect to coordinate-coordinate 
coupling, the free energy has two differences, namely, 
(i) the appearance of the extra factor G^ in the term 
AiFq{T,B), and (ii) the presence of the additional term 
A2Fo(T,B). It is interesting to see that in the absence 
of the magnetic field, the free energy has the same ex- 
pression irrespective of the form of coupling. 

With respect to the two schemes of coupling, the dif- 
ferences observed here in the form of the thermodynamic 
function, i.e., the free energy, add to the ones noted in the 
reduced dynamical description of the charged particle by 
means of the quantum Langevin equation [4:] . This work 
thus re-emphasizes the crucial role played by the form of 
the coupling with the heat bath in characterizing the be- 
havior of the charged particle. Investigations to obtain 
the different thermodynamic properties of the charged 
particle, e.g., its entropy and specific heat, by using Eq. 
([75)1 are in progress. 
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